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Abstract: The Anti-de Sitter-Reissner-Nordstrom (AdS-RN) black hole in the canon¬ 
ical ensemble undergoes a phase transition similar to the liquid-gas phase transition, i.e. 
the isocharges on the entropy-temperature plane develop an unstable branch when the 
charge is smaller than a critical value. It was later discovered that the isocharges on 
the entanglement entropy-temperature plane also exhibit the same van der Waals-like 
structure, for spherical entangling regions. In this paper, we present numerical results 
which sharpen this similarity between entanglement entropy and black hole entropy, by 
showing that both of these entropies obey Maxwell’s equal area law to an accuracy of 
around 1 %. Moreover, we checked this for a wide range of size of the spherical entan¬ 
gling region, and the equal area law holds independently of the size. We also checked 
the equal area law for AdS-RN in 4 and 5 dimensions, so the conclusion is not specific to 
a particular dimension. Finally, we repeated the same procedure for a similar, van der 
Waals-like transition of the dyonic black hole in AdS in a mixed ensemble (hxed electric 
potential and hxed magnetic charge), and showed that the equal area law is not valid in 
this case. Thus the equal area law for entanglement entropy seems to be specihc to the 
AdS-RN background. 
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1 Introduction 

Entanglement entropy (EE) appears to be a versatile tool that can be used to study a 
rich variety of physical phenomena. In particular, it can serve as a probe of the different 
phases of the theory [1, 2], ranging from the confining phase of large-A gauge theories 
[3] to topological phases in condensed matter systems [4], to tachyon condensation [5] 
and superconducting phase transitions [6]. 

Entanglement entropy has also emerged as a central component of the AdS/CFT 
correspondence. According to the Ryu-Takayanagi formula [7, 8], the entanglement 
entropy Sa between a boundary region A and its complement is computed (in static 
backgrounds) in an elegant geometric fashion as the area of a minimal surface. The 
striking similarity between the Ryu-Takayanagi formula and the Bekenstein-Hawking 
formula for black hole entropy suggests some deep connection between entanglement 
entropy and black hole entropy. It has even been suggested that the origin of black hole 
entropy is entanglement entropy [9-11]. 
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Motivated by the themes above, in this paper we track entanglement entropy across 
a family of van der Waals-like phase transitions of charged black holes in AdS. The hrst 
phase transition under study is the one of AdS-RN in the canonical ensemble in 4 and 5 
dimensions. This transition was first discovered in [12, 13]: the curves of constant charge 
on the temperature-entropy plane have an unstable portion when the charge is smaller 
than a critical value. Moreover, at critical charge, the unstable portion squeezes to an 
inflection point. It was subsequently pointed out in [14] that the same qualitative behav¬ 
ior can be observed if we study the isocharges in the entanglement enfropy-temperature 
plane (for spherical entangling regions), and that moreover the inflection point occurs at 
the same critical temperature as the black hole with the same critical exponent. These 
hndings were then generalized to a wider class of supergravity backgrounds in [15], where 
it was found that, in all cases, the isocharges on the entanglement enfropy-temperature 
plane mimick the qualitative behavior of the ones on the entropy-temperature plane. 

In this paper, we ask the question of how far we can push this similarity between 
the two kinds of entropies with respect to these phase transitions. In particular, we 
will investigate whether Maxwell’s equal area construction, which is known to hold in 
temperature-entropy plane, is also valid in the temperature-entanglement entropy plane. 
Maxwell’s construction in the context of black hole thermodynamics has generated some 
recent interest: this topic was studied in [16, 17] in the context of AdS-RN, but also in 
[18, 19] in other contexts. We will show that, for AdS-RN, the van der Waals behavior on 
the entanglement entropy-temperature plane indeed also obeys the equal area law, with 
the transition temperature obtained by minimizing the black hole free energy function. 
We also vary the size of the spherical entangling region over a wide range and checked 
that the conclusion is valid to an excellent accuracy. Moreover, we checked this for the 
4-dimensional and 5-dimensional black holes, so that the equal area law is not specific 
to a particular dimension. We also repeated the same procedure for another van der 
Waals-like phase transition: that of the dyonic black hole in the hxed electric potential 
and fixed magnetic charge ensemble, and our numerical results indicate that the equal 
area law does not hold in this case. Thus, the equal area law seems to be a specihcity 
of the AdS-RN background. 

The rest of the paper is organized as follows: In Section 2, we review the phase struc¬ 
ture of AdS-RN in 3-1-1 dimensions in the canonical ensemble and discuss Maxwell’s equal 
area law in the entropy-temperature plane. In Section 3, we then turn to the numerical 
computation of holographic entanglement entropy, and present numerical evidence that 
the equal area construction also holds on the entanglement entropy-temperature plane. 
In Section 4, we repeat for the dyonic black hole in 3-1-1 dimensions and show that the 
equal area construction is not valid for this background. Next, in Section 5, we check 
that the equal area law also holds for AdS-RN in 4-1-1 dimensions. Finally, in Section 6, 
we summarize our main hndings and discuss a few possible future work. 


2 


2 Review of 4d AdS-RN in the canonical ensemble 


In this section, we survey the phase structure of the 4-dimensional AdS-RN black hole 
in the hxed charge ensemble, leading up to the van der Waals behavior in the entropy- 
temperature plane (i.e. there exists a family of hrst order transition ending with a second 
order one) and Maxwell’s construction The Einstein-Maxwell action in 4 dimensions 
reads: 

/ = I d^xV^iR - 2A - F^). (2.1) 

The AdS-RN solution together with the gauge held is given by: 

ds^ = -f{r)de + -7^ + r\de^ + sin^ 0d02), (2.2) 

f{r) 


fir) 


1 


2M Q2 ^2 

-L -5-1- 

r 


A = Q ( - ] dt. 

r J 


(2.3) 

(2.4) 


where M is the mass, Q is the electric charge and L is the AdS lengthscale. The additive 
constant in At was chosen to be — so that the norm of the vector potential A^ is regular 
at the horizon. The black hole temperature and entropy are: 


3< + L\rl - Q^) 
47rL^rf 


(2.5) 


S' = vrr^ , (2.6) 

where r_|_ is the horizon (the largest root of /(r_|_) = 0). From (2.5) and (2.6), we can 
easily eliminate the parameter r_|_ to obtain the function T(S, Q): 


T{S,Q) 


1 

47r 




-Q^ 


7r3/2\ 

S3/2 j ■ 


(2.7) 


From the function T{S, Q) above, one can plot the isocharges on the T — S plane. The 
plot is presented on the right panel of Figure 1. As can be seen from the plot, the curve 
is monotonic for sufficiently large Q. As Q decreases, the curve has an inflection point 
when Q reaches a threshold value Qc- One can solve for the position of the inflection 
point by: 


dT\ _ rd^T\ 


( 2 . 8 ) 


^The distinction between “first order” and “second order” here refers to the slope of the free energy 
plot versus the temperature. Upon a closer look, the nature of these phase transitions is more subtle: 
see for example [20] where the phase transition was studied using Ehrenfest equations. 
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We find the critical entropy critical charge Qc and critical temperatnre to be: 


= -L^ , 

6 

(2.9) 

L 

Qc = g ’ 

(2.10) 

[2 1 

(2.11) 



Finally, when Q < the cnrve becomes oscillatory, and there is a small portion with 
negative heat capacity: 

r(|)^<0, (2.12) 

Like in the case of the liqnid-gas transition, this portion is believed to be thermodynam¬ 
ically nnstable, and shonld be replaced by an isotherm T = T* according to Maxwell’s 
prescription. The valne of T* can be obtained in two different (but equivalent) ways: by 
the equal area condition, or from the Helmholtz free energy. The hrst method, the equal 
area condition, states that T* is the unique temperature which divides the oscillatory 
part of the curve T{S) into two regions with equal area. We will find T* from the second 
method, i.e. using the Helmholtz free energy, and check numerically that it is equivalent 
to the first. The Helmholtz free energy can be found from the on-shell action 

1 / , , 3Q‘^L^\ 


We present in the left panel of Figure 1 the plot of F versus T. For Q < Qc, we 
observe the swallowtail behavior familiar from catastrophe theory, and the transition 
temperature T* is the horizontal coordinate of the junction between the two stable 
branches. Numerically, using Q = 1.5 and L = 10, we found T* 0.02663. We can now 
proceed to verify Maxwell’s equal area law, which is the following statement: 


-52 


T{S, Q)dS - %{S 2 - S^) = T,(^3 - ^ 2 ) 


-53 


'Si 


'S 2 


T{S,Q)dS 


(2.14) 


where Si, S 2 and S 3 are the smallest, intermediate, and largest roots of the equation 
T* = T{S,Q). Numerically, we found that both sides of (2.14) evaluate to around 
0.00761, thus confirming the validity of Maxwell’s construction. 

^The Helmholtz free energy of AdS-RN is usually measured with respect to an extremal background 
with the same electric charge [12]. For our purposes, however, this background subtraction only shifts 
the plot of F versus T in the vertical direction and does not affect the transition temperature T*. 
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Figure 1: Plot of the temperature versus the free euergy (left pauel) aud plot of the 
temperature versus the eutropy (right pauel). For both pauels, the values of the charge 
are Q = 1.5 (greeu), Q = Qc = y aud Q = 1.8 (orauge), aud we use L = 10. The 
trausitiou temperature for the greeu curve is T* = 0.02663. 

3 Maxwell construction for entanglement entropy 

lu this sectiou, we iuvestigate the phase structure of the EE-temperature plaue iustead 
of the eutropy-temperature plaue. As discussed iu the iutroductiou, the work of [14] 
demoustrates that, at least for AdS-RN iu 4d, the isocharges ou the EE-temperature 
plaue mimick the qualitative behavior of the eutropy-temperature plaue. However, that 
paper did uot make the couuectiou with Maxwell’s coustructiou. Iu this sectiou, we 
will show that, iudeed, eutauglemeut eutropy also respects the Maxwell coustructiou, 
thereby streugtheuiug the couclusiou of [14]. First, let us briefly review a few held 
theory geueralities about EE. 

3.1 Brief review of holographic EE 

Suppose we have a quautum held theory described by a deusity matrix p, aud let A be 
some regiou of a Cauchy surface of spacetime. The eutauglemeut eutropy betweeu A 
aud its complemeut A”’ is dehued to be: 

Sa = -Tr^(pAlogp^), (3.1) 

where pA is the reduced deusity matrix of A\ pA = Tr^c(p). As meutioued iu the 
iutroductiou, eutauglemeut eutropy is computed holographically by the Ryu-Takayauagi 
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recipe 

^ _ Area(r^) 

where is a codimension-2 minimal surface with boundary condition c^Fyi = dA^ and 
Gat is the gravitational Newton’s constant. 

A few remarks are in order concerning EE at hnite temperature and in hnite volume. 
At nonzero temperature, entanglement entropy no longer has the nice properties that 
it does at zero temperature, for example the area law (see, for example, [21]). This is 
due to the fact that, at hnite temperature, entanglement entropy is “contaminated” by 
a thermal component which scales as the volume of the entangling region rather than 
its area. 

Moreover, when the bulk is topologically nontrivial (as is the case for AdS-RN with 
spherical horizon), the Ryu-Takayanagi formula is rehned by an additional topological 
constraint: only surfaces which are homological to the entangling region on the boundary 
are considered in the minimization problem [22]. This constraint, which ensures that the 
Araki-Lieb inequality is satished, implies that, for sufficiently large entangling region, the 
minimal surface is disconnected and includes the horizon itself as a connected component. 
It might appear curious that the homology constraint means the entanglement entropy 
of A is not equal to that of the complement of A, but this is expected when the system 
is in a mixed state. By avoiding large entangling regions, we will also avoid having to 
deal with this “phase transition” between connected and disconnected minimal surfaces 

4 


3.2 Maxwell’s equal area law 

We will take the region A to be a spherical cap on the boundary delimited hy 9 < 9 q. 
From the remarks in the previous subsection, we will pick four values of 9q which are 
smaller than 9q = 0.1,0.15,0.2 and 0.4. The minimal surface can be parametrized 
by the function r{9), which is independent of the coordinate (j) by rotational symmetry. 
The area functional is: 


A 


2t\ 


r-00 


r sin^ 



-I- r‘^d9 , 


(3.3) 


where r' = ^. The function r[9') is the obtained by solving the Euler-Lagrange equation: 

dL d fdL\ 
dr d9 \5r'/ 

^The Ryu-Takayanagi formula only applies to static backgrounds and when the bulk theory is Ein¬ 
stein gravity. 

^Similar phase transitions are also observed in infinite volume, when the entangling region is com¬ 
posed of multiple strips [23]. 


(3.4) 
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with the boundary conditions r( 6 'o) —)■ oo and r'( 0 ) = 0 (i.e. the minimal surface is 
regular at the center 6 ^ = 0 ). 

Also, since EE is UV-divergent, it has to be regularized. We will do it by subtracting the 
area of the minimal surface in pure AdS whose boundary is also 9 = 9q. In other words, 
we hrst integrate the area functional to some cutoff 6 *c < 6 *o. Then we set M = Q = 0 
to obtain AdS in global coordinates: 

ds^ = -(1 + -)de + --^ + T^dnl. (3.5) 

^ 1 + X2 

The minimal surface which goes to 9 = 9 q on the boundary is given by 

( 2 \ “h2 

(S) - ■) « 

We can easily integrate the area functional of this minimal surface up to 9^. Then we 
subtract this quantity from the black hole one to obtain the renormalized entanglement 
entropy, which we will denote by If we could compute the area of the minimal 

surface analytically, then we have to take the limit 9c —)■ 6^0 after the subtraction. But 
since we are doing a purely numerical computation, all we can do is choose 9c close to 

9o- 

In Figure 2, we present the plots of the isocharges on the AS'^ — T plane for the 4 values 
of 6 * 0 . In the case 6 *o = 0.1, we present 3 isocharges corresponding to a charge below the 
critical one, equal to the critical one, and above the critical one. As can be seen on this 
plot, the van der Waals phase structure noted in [14] is observed. For the 3 other values 
of 6 * 0 , only a subcritical isocharge is presented since this is enough to verify the equal 
area law. We have also drawn the transition isotherm in dashed green taken from the 
free energy function. Finally, we proceed to check whether Maxwell’s equal area holds. 
We will refer to the two areas to be calculated as Ai and A 2 : 

Ai= / T{ASA,Q)dASA-n{AS^A^ - (3.7) 

JAsy 

pAsy 

A 2 = ryASf - A^f) - / T(AAa, Q)dASA (3.8) 

jAsy 

where AS^2\ AS^^'^ and AS^^'^ are the smallest, intermediate, and largest roots of the 
equation T* = T{ASa, Q)- Maxwell’s equal area would amount to the statement: 

Ai = A 2 (3.9) 

We tabulate in table 1 the values of Ai and A 2 for each choice of 6 * 0 , as well as the 
choice of the cutoff 6 *c for each and the relative error between Ai and A 2 (taken to 
be the difference between Ai and A 2 divided by their average). As can be seen from 
the small relative errors in table 1 , it is safe to claim that the equal area law holds for 
AdS-RN in 4 dimensions. 
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Figure 2: Plot of isocharges on the T — ASa plane, for 9q = 0.1 (top left), 0.15 (top 
right), 0.2 (bottom left) and 0.4 (bottom right). For all panels: the values of the charge 
are Q = 1.5 (green), Q = Qc = ^ (red) and Q = 1.8 (orange). The transition isotherm 
is obtained from the free energy (left panel of Figure 1). For the green curves, we also 
show the data points which were used to create the interpolation. 


4 A counterexample: dyonic AdS-RN 


In this section, we repeat the procedure above for another charged black hole in AdS: 
the dyonic solution in 4d, which describes a black hole in AdS carrying both an electric 
charge and a magnetic charge. This black hole is also a solution to the Einstein-Maxwell 

















^0 

Oc 

Ai 

A2 

Relative error 

0.1 

0.099 

1.59197 X 10“® 

1.56443 X 10-® 

1.74 % 

0.15 

0.149 

5.40162 X 10-® 

5.30496 X 10-® 

1.8 % 

0.2 

0.199 

0.0000127702 

0.0000126147 

1.22 % 

0.4 

0.399 

0.000104962 

0.000103373 

1.52 % 


Table 1: Comparison of Ai and A 2 for the AdS-RN black hole in 3+1 dimensions. 


action (2.1). The line element together with the gauge held are given by [24-26]: 


ds^ 


-f{r)df + 


dr"^ 

f{r) 


+ r‘^{d6‘^ + sin^ 6d(p‘^) , 


(4.1) 


with 


A 



^ dt -h P cos 9d(j) , 


fir) 


2M 

1-h 

r 


Q^ + P^ 



Here P is the magnetic charge. 


(4.2) 

(4.3) 


4.1 The van der Waals transition of the “mixed” ensemble 

As for the AdS-RN solution, the phase structure depends on which statistical ensemble 
one chooses. One can, for example, study the ensemble with both the electric and 
magnetic charges hxed. This is however not very interesting, since the results in this case 
can be trivially obtained from the results in section 2 with the replacement —)■ Q^+P^. 
Thus, we will instead work in the hxed <h, hxed P ensemble. It was observed in [26] 
that the phase structure in this “mixed” ensemble exhibits the van der Waals transition. 
The asymptotic value of the electric potential <h is related to the electric charge Q by: 

$ = — . (4.4) 

r+ 


The temperature and entropy are given by: 


T 



(4.5) 


S = vrr^ . (4.6) 

From which we hnd the function T(5, <h,P). Also, since the electric charge is now 
allowed to vary, we should compute the Gibbs potential: 


G = M -TS 


(4.7) 
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Figure 3: Plot of the temperature versus the Gibbs potential (left) and the temperature 
versus the entropy (right) for the dyonic black hole. In both panels, the values chosen 
for the magnetic charge are P = 0.95 (green), P = Pc = 1.0667 (red) and P = 1.3 
(orange). We also set $ = 0.6 and L = 10. The transition temperature for the green 
curve is T* = 0.02135. 


instead of the Helmholtz free energy. Computation gives: 


G 


3P2 

4r+ 


W 

4 


ti - 



(4.8) 


If we keep $ hxed, and plot the curves of constant P on the S — T plane, we can observe 
a van der Waals-like transition when P is smaller than a critical value Pc- The plot of 
these curves is presented on the right panel of Figure 3. The critical magnetic charge 
Pc-, entropy Sc and temperature Tc can be found to be: 



(4.9) 




(4.10) 

SLn 

(4.11) 


We plot the Gibbs potential on the left panel of Figure 3. From this plot, we obtained 
the transition temperature T* = 0.02135 with P = 0.95, $ = 0.6 and L = 10. We then 
checked numerically that T* is indeed the temperature which obeys Maxwell’s equal area 
law: numerical integration of both sides of (2.14) yields a value of around 0.004658. 
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Figure 4; Plots of curves of constant magnetic charge on the T — AS a plane for the 
dyonic black hole in the “mixed” ensemble, with do = 0.1 (left) and 0.15 (right). For both 
panels, the chosen values of the magnetic charge are P = 0.95 (green), P = = 1.0667 

(red) and P = 1.3 (orange). We set <F = 0.6 and L = 10. In both plots, the dashed 
isotherm is obtained from the Gibbs potential (left panel of Figure 3). For the green 
curves, we also show the data points which were used to create the interpolation. 


Oo 

dc 

Ai 

A2 

Relative error 

0.1 

0.15 

0.0995 

0.1499 

1.60959 X 10-® 
5.467 X 10“® 

1.25275 X 10-6 
4.2505 X 10-6 

24.93 % 

28.62 % 


Table 2: Comparison of Ai and A 2 for the dyonic black hole. 


4.2 Maxwell’s equal area law for entanglement entropy 


Next, we turn to the computation of entanglement entropy, and we will show that 
this time the equal area law is not valid for entanglement entropy, unlike black hole 
entropy. We will consider two values of Oq: 0.1 and 0.15. The plots of temperature 
versus entanglement entropy at fixed P are presented in Figure 4 (with $ = 0.6 and 
L = 10). Again, in table 2 we tabulate the values of the two areas Ai and A 2 for each 
choice of do, as well as the choice of dc and the relative error between Ai and A 2 . As 
can be seen from the large relative errors in table 2, the equal area law does not seem 
to hold for the dyonic black hole in the mixed ensemble. 
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5 Another example: 5d AdS-RN 

In this section, we will show that the eqnal area law for entanglement entropy also works 
for AdS-RN in 4+1 dimensions. 


5.1 The van der Waals transition in the canonical ensemble 

The AdS-RN solution in arbitrary dimension is given in [12], For 4+1 dimension, the 
metric is: 

ds^ = —fir)dt^ + + sin^ '0(d6'^ + sin^ ddcf)^)) (5.1) 

fyr) 

where are hyperspherical coordinates on the 3-sphere, with ip and 6 ranging 

from 0 to TT, and (j) ranging from 0 to 27r, and 


■’ ^ ’ 37rr2 ^ 37r2r4 ^ 

The black hole temperature and entropy are given in terms of the horizon r+ by: 

^ r+ 1 2Q^ 

7rL2 27rr_|_ 


(5.2) 


(5.3) 


S = yri (6.4) 

From which we obtain the function T(S', Q). Following the same steps as the 4-dimensional 
case, we plot the isocharges on the T — S plane on the right panel of Figure 5. We find 
the inflection point to be at: 



(5.5) 


(5.6) 

~ _ 4V3 

^ 57rL 

(6.7) 


The free energy in the fixed charge ensemble is now: 


F 


71 



rt + 


20Q^\ 

37r2r+ / 


(5.8) 


We plot the free energy versus the temperature on the left panel of Figure 5. The 
phase structure again presents a van der Waals transition, like in 4d. The transition 
temperature was found to be T* = 0.149086. We can now check numerically the equal 
area law: evaluating both sides of (2.14) numerically yields around 0.01711. 
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Figure 5: Left panel: Plot of the free energy versus the temperature for AdS-RN in 
4+1 dimensions, with L = 3 and Q = 1. Right panel; Plot of isocharges on the T — S 
plane also with L = 3 and Q = 1. The transition temperature for the green curve is 
T* = 0.149086. 

5.2 Maxwell equal area law for entanglement entropy 

Next, we turn to computing the entanglement entropy of a spherical region. The 5d 
analog of the disk-shaped region 6 < 6o in 4d is the spherical region 'ip < 'ipQ for some 
constant -00 ^ [Oj^r], We will choose 4 different values of 'ipQ- 0.05, 0.1, 0.15, and 0.8. 
The minimal surface can be parametrized by r('^), where this function does not depend 
on 0 or 0 by rotational symmetry. The area functional is given by: 



(5.9) 


Proceeding as for the 4d case, we will solve the equation of motion numerically with the 
boundary condition that the minimal surface is regular at the center and goes to 'ipQ at 
the boundary. To regularize entanglement entropy, we again integrate to some cutoff 'ipc 
and subtract the pure AdS minimal surface which goes to the same 'ipo at the boundary. 
Such a surface is analytically given by: 



(5.10) 


We present in Figure 6 plots of isocharges for each choice of "00 • As previously, we 


tabulate in table 3 the values of Ai and A 2 for each case, and the relative error between 
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Figure 6: Plot of isocharges on the T — ASa plane for the 4+1 dimensional AdS-RN 
solution, and -00 = 0-05 (top left), 0.1 (top right), 0.15 (bottom left) and 0.8 (bottom 
right). The values of the charge chosen are Q = 1.75 (green), Q = Qc = (red) and 
Q = 3 (orange). We have set L = 3. The transition isotherm is obtained from the free 
energy (left panel of Figure 5). We also show the data points which were used to create 
the interpolation. 


Ai and A 2 and the choice of the cutoff value 'ipc- As in the 4d AdS-RN case, the relative 
errors are of the order of 1%, and therefore we can safely claim that the equal area law 
works. 
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"00 

i’c 

Ai 

A 2 

Relative error 

0.05 

0.0499 

1.12124 X 10“^ 

1.10842 X 10"^ 

1.15 % 

0.1 

0.0999 

4.05257 X 10“^ 

3.98178 X 10-^ 

1.76 % 

0.15 

0.1499 

2.06799 X 10“® 

2.05255 X 10“® 

0.75 % 

0.8 

0.799 

0.00176189 

0.00171273 

2.87 % 


Table 3: Comparison of Ai and A 2 for the AdS-RN black hole in 4+1 dimensions. 


6 Conclusion 

The laws of black hole thermodynamics have provided us with a robust understanding 
of black holes as thermodynamical systems, though the nature of the microscopic con¬ 
stituents of these systems remains by and large mysterious. The gauge-gravity duality 
seems to suggest, on very general ground, that gravity could emerge from the dynamics 
of a strongly coupled large-A quantum held theory. In view of this, it is an important 
thing to do to see whether observables in such a quantum held theory can mimick the 
behavior of gravitational systems. 

In this paper, we have presented compelling numerical evidence that, while a large 
class of charged backgrounds in AdS exihibits van der Waals-like phase transitions (in 
an appropriate statistical ensemble), the transition of the AdS-RN background is special 
in that its holographic EE obeys the Maxwell’s equal area law (for spherical entangling 
regions). In this respect, we have improved upon previous studies such as [14]. More¬ 
over, the equal area law for entanglement entropy seems valid regardless of the size of 
the entangling region and the dimensionality of the background. That said, we have 
only considered spherically symmetric entangling regions, and future work to generalize 
the conclusions in this paper to entangling regions of other shape could yield further 
interesting insights. 

Finally, we note that the van der Waals transition of AdS-RN has been observed in 
the context of the extended black hole thermodynamics, where the cosmological constant 
is identihed as the pressure variable and its thermodynamic conjugate as the volume 
variable [27]. In this framework, one can analyze the isotherms on the P — V plane, 
and they turn out to be remarkably similar to the van der Waals gas. One can wonder 
about the relationship between the van der Waals transition in the P — V plane and the 
one in the T — S plane. It was pointed in [16] that the two are related to each other 
by a duality similar to T-duality of string theory. Moreover, it was noted by [19] that 
Maxwell’s construction in the P — V plane can be surprisingly subtle (the construction 
does not work if the volume is replaced by the specific volume, unlike the usual van der 
Waals gas). We note that the holographic interpretation of the P — V plane is not very 
well-understood. However, the van der Waals transition on the P — V plane seems to 
be connected with the renormalization group flow (see [15, 28]). 
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Finally, it would be interesting to extend our investigations to other charged black 
holes in AdS, for example (among others) to more complicated conhgurations with mag¬ 
netic held. For instance, there has been lots of recent efforts to construct magnetic stars 
in AdS (see, among others, [29-31]) which are the gravitational duals to condensed mat¬ 
ter systems. In particular, entanglement entropy has been studied in such a background 
in [32]. Even though the equal area law most probably will fail in such backgrounds, 
the study of other backgrounds may shed light on why it works in the hrst place for 
AdS-RN. 
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